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;h ' Abstract 



We consider the production of a Higgs boson via gluon-fusion and its decay into 
two photons. We compute the NLO virtual QCD corrections to these processes in a 
general framework in which the coupling of the Higgs boson to the external particles 
is mediated by a colored fermion and a colored scalar. We present compact ana- 
lytic results for these two-loop corrections that are expressed in terms of Harmonic 
Polylogarithms. The expansion of these corrections in the low and high Higgs mass 
regimes, as well as the expression of the new Master Integrals which appear in the 
reduction of the two-loop amplitudes, are also provided. For the fermionic contribu- 
tion, we provide an independent check of the results already present in the literature 
concerning the Higgs boson and the production and decay of a pseudoscalar particle. 
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1 Introduction 



The Higgs searches program at the TEVATRON and at the LHC requires from the theo- 
retical side the highest possible level of accuracy in the prediction of the production cross- 
sections and of all the decay channels. Over the years a lot of effort has been devoted to 
the study of the QCD, and also EW, corrections to the various production mechanisms and 
decays in the Standard Model and beyond (for a recent review see Ref. [1]). 

The gluon-fusion process gg H + X [2J is the dominant production mechanism. Its 
present knowledge includes the NLO [H HI [5] and NNLO QCD corrections [6] and the two- 
loop EW corrections [HlHlin]- The QCD corrections to Higgs production at finite transverse 
momentum have also been discussed [ID]. While the NLO QCD corrections and the two- 
loop EW light fermion contribution are known completely, namely for arbitrary value of the 
Higgs mass and of the other relevant particles in the loops, the NNLO QCD corrections are 
only known in the heavy top limit while the result for the two-loop EW top contribution is 
valid only for intermediate Higgs mass, i.e. ttIh < '2,mw 

The Higgs decay if — 77 [JJj| is, for light values of the boson mass, a very promising 
channel. It has been studied in great detail including the NLO QCD [121 [13] and the two- 
loop EW corrections [El [HI [15], [16] . The NLO QCD corrections are now known in a closed 
analytic form [T3l [5], while for the EW corrections their knowledge is similar to that of the 
gluon fusion process. 

Given the importance of the Higgs physics program, it is highly desirable to have the 
radiative corrections to the various reactions expressed in analytic form that can be eas- 
ily implemented in computer codes. With respect to this, it should be recalled that the 
complete result concerning the NLO QCD corrections to the gluon fusion process has been 
reported in Ref. 0] via a rather lengthy formula expressed in terms of a one-dimensional in- 
tegral representation. Actually the calculation of the two-loop light-fermion EW corrections 
to the Higgs production and decay ^ has shown that corrections of this kind can be calcu- 
lated analytically, expressing the results in terms of Harmonic Polylogarithms (HPL) [IT], 
a generalization of Nielsen's polylogarithms, and an extension of the HPL, the so-called 
Generalized Harmonic Polylogarithms (GHPL) [18]. The idea lying behind the introduc- 
tion of (G)HPLs is to express a given integral coming from the calculation of a Feynman 
diagram in a unique and non-redundant way as a linear combination of a minimal set of 
independent transcendental functions. These functions are expressed as repeated integra- 
tions over a starting set of basis functions and this set depends strongly on the problem 
one has to solve, being connected directly to the threshold structure of the diagrams under 
consideration. 

An inspection of the threshold structure of the NLO QCD corrections to the gluon- 
fusion process and to — 77 decay shows that these corrections can be fully expressed 
in terms of the original set of HPLs introduced in [17] . A FORTRAN program [19] and a 
Mathematica package [20] that efficiently evaluate these functions are available. 

The aim of this paper is to provide analytic expressions, in terms of HPLs, for the NLO 
QCD corrections to the Higgs production cross section via gluon fusion, i.e. gg — if, in a 
general form that can be applied both to the SM and to models beyond it, and, moreover, 
to provide an independent check for the formulas already present in the literature. The 
production mechanism is assumed to be mediated by colored fermion and scalar particles. 
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As a byproduct we also present the NLO QCD corrections to the Higgs decay into two 
photons, i.e. H 77. A similar project has been carried out in Ref. [5j. There, the authors 
started from the result of Ref. [4]^ expressed as a one-dimensional integral representation. 
Expanding this result in a power series, employing the theorem that two analytic functions 
are the same if their Taylor series are the same, they were able to rewrite it in terms of 
HPLs. In our case we explicitly compute all the relevant Feynman diagrams, expressing 
the result in term of HPLs. The calculational techniques we employed are the Laporta 
algorithm [21] for the reduction to Master Integrals (Mis) and the differential equation 
method [22] for their calculation (the calculation is implemented in FORM [23] codes). To 
complete an independent check of the results presented in Ref. |4] we also computed the 
NLO QCD corrections to the pseudoscalar production and decay, i.e. gg ^ A, A ^ 77. 

The paper is organized as follows: in Section [2] we discuss the QCD corrections to the 
decay width H 77. Section [3] is devoted to the study of the Higgs production via the 
gluon fusion mechanism. The following section contains the analytic expressions for the 
virtual QCD corrections to the fermionic contribution in gg A, A 77. Finally we 
present our conclusions. We include also two Appendices. In the first one we collect the 
expansions of the relevant functions in the two regimes: for Higgs mass much lighter than 
the particles mediating the Higgs interaction with the vector bosons and in the opposite 
case. In the second Appendix we collect the Mis not already present in the literature, that 
enters the calculation of the NLO QCD corrections. 



2 The ^ 77 Decay Width 

We begin by considering the decay width H 77. Being the Higgs boson electrically 
neutral its coupling to the photon is mediated at the loop-level by charged particles. For 
the latters we assume a vector boson neutral under SU{Nc), a fermion and a scalar particle 
in a generic R1/2, Ro SU{Nc) representation, respectively, whose coupling's strengths to the 
Higgs are: 

HVV = gXim^, HFF = gXy2^, HSS = gXo—, (1) 

where g is the SU{2) coupling, m^^ is the W mass, mi/2 is the fermion mass, A is a generic 
coupling with the dimension of mass and Aj are numerical coefficient^. 
The partial decay width for the reaction if — 77 can be written as: 

m - 77) = ^^^^ i-^i' , (2) 

where the function can be organized with respect to the lowest order term and its QCD 
corrections as: 

T = \iQl Ni + Ai/2 Q?/2 + Ao Ql Nq— , (3) 



^In Ref.jl] the QCD corrections were considered only for the fermion contribution. 
^The SM is recovered with Ai = A1/2 = 1, Aq = 0, Ai'c = 3 and R1/2 = 3. 
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where tuq is the mass of the scalar particle, while Qi and iVj, i = 0, 1/2, 1, are the electric 
charges and the representation numbers under SU{Nc) of the scalar, fermion and vector 
boson particles, respectively. 
Writing: 

we have at the one-loop level 

^f'^ = 2(l + 6yi)-12yi(l-2i/i)i7(0,0,xi), (5) 

•^V? = -%i/2 [2 - (1 - %i/2) H{0, 0, X1/2)] , (6) 

J-f^ = 4yo[l + 2 2/oi/(0,0,xo)] . (7) 

In Eqs.dMZD 

_ _ VI - %i - 1 /ox 

^• = ^' "•= vi-% + i - 

with mi the mass of the vector particle and, employing the standard notation for the HPLs, 
if(0, 0, z) labels a HPL of weight 2 that results to bcl 

H{Q,0,z) = ^\og'{z) . (9) 



The QCD corrections to the lowest order result can be written as 

as 

7T 



-^a = v E (10) 

i=(0,l/2) 

where C (R) is the Casimir factor of the Ri representation (in particular, for the fundamental 
and the adjoint representations of SU{Nc) we have Cp = {N^ — l)/{2Nc) and Ca = Nc, 
respectively). We consider first the fermion contribution (the relevant Feynman diagrams 
are shown in Fig. [T] (a)-(d)). 

The expression for J-'[^^2 depends upon the renormalized mass parameter employed. In 
the case of MS quark masses we have 

= + 1" (^) . (11) 



where is the 't Hooft mass and 

W2«,a). X _ 36x Ax (1 - 143; + x^) 4x{l + x 
y"") - (a;_i)2 (a; -1)4 (x - 1)3 



8x (1 + 9a; + x'^) ^^.^ ^ , 2x (3 + 25a; - 7x'^ + 3x^^ 



H{0, 0, x) + — ^ -. H{0, 0, 0, x) 



(a; -1)4 ^ ' ' ^ (x-l)5 

+ fa _ 1)1 IC^-ff (0. ^) + * H(.o, -1, 0, 1) - (0, 1, 0, x)] 

^ 4x(5 -6x + 5x^) ^^^^^^^^^^ _ 8.(1 +^:^_+^^^ + ^') ^,(,) ^ (12) 



•^AU the analytic continuations are obtained with the replacement — 
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U) (k) {1} {m) (n) 



Figure 1: The Feynman diagrams for the decay process H,A^ 77. Diagrams (a)-(d) have 
a fermion, labeled by "f", running in the loop, while diagrams (e)-(n) have a scalar, labeled 
by "s". 



with 



9.2 ^ ^ ^ 1 . 7 



n^ix) = -C2' + 2C3i/(0,x) + C2i^(0,0,x) + -i/(0,0,0,0,x) + -i/(0,l,0,0,x) 

-2 H{0, -1, 0, 0, x) + 4 H{0, 0, -1, 0, x) - H{0, 0, 1, 0, x) . (14) 

In Eqs.f ll2|[T4|) = C(^) the Riemann's zeta functions. 

The expression for ^-"[^^2 in case the one-loop result is given in terms of on-shell fermion 
masses is given instead by: 

:Fi^^ = :F[^r\xy,) + l:F[^^\xy,). (is) 

Eq. ffTSl) is in agreement with the results presented in [131 IS]- 

(20 

We now present the scalar contribution, jFg , assuming that both the mass of the scalar, 
mo, and the coupling A are renormalized in the MS scheme (the relevant Feynman diagrams 
are shown in Fig. [T] (e)-(n)). We find 

= "^^(^0) + (^f '^^(xo) + '^^xo)) In , (16) 



where 



Ux 24x2 a; (3 - 8x + 3^^) 34^^ 

(x-i)2 " (x-i)4^'+ (x-i)3(x + i) ^ ''''^J!^:^'' 
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_ -- [C2H{0, x) + 4H{0, -1, 0, x) - H{0, 1, 0, x) + H{1, 0, 0, x)] 
[X ij 

2x^(5 -11a;) , , 16x^(1 + ^2) 

= 7 ..3. ^,, g(0,x)-- — if(0,0,x), (18) 

(x — l)'^(x + 1) (x — Ij^ 

^r)(x) = -^^r- (19) 

(21) 

We provide also JFq assuming that the mass of the scalar is renormalized on-shell while 
the A coupling is still given as an MS one. It reads 



= + I :Fr\^o) + -^r'^^(xo) in . 



(20) 



3 Virtual Corrections to gg ^ H Production Mecha- 
nism 

In this section we present the analytic expressions for the virtual two-loop QCD corrections 
for Higgs boson production via the gluon fusion mechanism. Being the Higgs boson neutral 
under SU{Nc), its coupling to the gluons is mediated by a loop of colored particles. As in 
Section [2], we consider a fermion and a scalar particle, that run in the loops. The Feynman 
diagrams relevant for the NLO corrections to the production cross section are shown in 
Fig.H 

The hadronic cross section can be written as: 

a{hi + h2^H + X) = / c?a;i(ix2 /a,/ii(xi,/i|) /b,?,2(x2,/i|) X 

„^ Jo 



X dz 5 (z Ciabiz) , (21) 

Jo V 2:1X2/ 

where th = w^l/s, fJ^F is the factorization scale, fa,hi{x, fijp), the parton density of the 
colliding hadron hi for the parton of type a, (a = g, q, q) and dah the cross section for the 
partonic subprocess ah ^ H + X a,t the center-of-mass energy s = X1X2S = m\/ z. The 
latter can be written as: 

a,fe(2) = a(°)zG,fe(z), (22) 

where 



128^2 



TT 



— n 1 /o V U / 



2 

(23) 



1=0,1/2 

is the Born-level contribution with Qf^^^ = J^^^^^ and T{Ri) are the matrix normalization 
factors of the Ri representation {T{R) = 1/2 for the fundamental representation of SU{Nc), 
T{R) = Nc for the adjoint one). 
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Figure 2: The Feynman diagrams for the production mechanism gg —>■ H,A. Diagrams 
(a)-(h) have a fermion, labeled by "f", running in the loop, while diagrams (i)-(x) have a 
scalar, labeled by "s". 



Up to NLO contributions we can write 



TT 



with 



= 5{1-Z)5ag5bg 

= S{l-z) 



+ 



(24) 

(25) 
(26) 



The dots in Eg. (126!) represent the contribution from the real emission that we have not 

written as well as the other NLO factors G'S(z), G^qq\z), (2). In Eq.([26]) /5o = (11 C^- 

4:nfT{Rf) — nsT[Rs)/6 with Uf {ug) the number of active fermion (scalar) flavor in the 

representation Rf {Rs)- 
(21) 

The function can be cast in the following form: 



6 



,(2/) 



rrir 



T{R,)(c{R,)gl'''''''\x,) + CAg.. 



,{21,Ca), 



Xi 



A' 



2 \ l-2j 



(10 



h.c. 



(27) 



with Qf'^'^^^ = J-'f'^\ The infrared regularized functions Qf'^''^'^\ after subtraction of the 
infrared poles, are found to be: 



{21,Ca). 
/2 ' 



Ax 



X] 



{x - ly 



3 + £(l±^^i/(0.0.0..)-i(l±i#H,W 



+ C3-i^(l,0,0,a;) 



(x-l) 



Ax 



{x - ly 



X] 



(28) 
(29) 



with 
'H2{x) 



C2' + 2C3 + ^ i^(0, a;) + 3C3i^(l, x) + (2^^ (1, 0, x) + ^ (1 + 2(2) H{Q, 0, x 



-2/7(1,0,0,3;) + iJ(0, 0,-1,0, x) + -if(0, 0, 0, 0, x) + 2 if(l, 0, -1, 0, x) 
-if(l,0,0,0,a;) . 



(30) 



The analytic expression of Q^'jl is in agreement with that reported in Ref.[5] based on the 
results presented in Ref . 0] . 



4 Pseudoscalar Higgs: A ^ 77 and gg ^ A 

To complete an independent check of the results of Ref. [1] in this section we consider the 
virtual NLO QCD corrections to the decay width of a pseudo-scalar particle A in two 
photons, T{A —> 77), and to its production cross section via gluon fusion, a{gg A). 

As in Ref. PJ we assume the interaction of the A particle with gluons mediated only 
by the top quark. Because the NLO QCD corrections to these two processes are calcu- 
lated in Dimensional Regularization, a prescription for the 75 matrix is needed. We use 
the same prescription of Ref. [5, i.e. 't Hooft-Veltman one [23], that, as it is well known, 
breaks manifestly Ward Identities. The latters need to be restored explicitly, with a finite 
renormalization. If Zj^jj and Z^^jj are the renormalization constants of the vertex Hff, 
scalar Higgs-fermion-antifermion, and Aff, pseudo-scalar Higgs-fermion-antifermion, re- 
spectively, the contribution of the finite renormalization can be found using [11 [25]: 

^Aff = ^Hff +'^Cf— ■ (31) 
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4.1 Decay Width ^ ^ 77 

In analogy with Eq.(j2]), we write: 



128y27r3 



r{A ^ 77) 



1^1 



(32) 



Assuming the strength of the couphng of the pseudoscalar to the top quark equal to Att = 
gritiTLt/ {2m^), with nit the top mass and rjt a numerical coefficient, the function S can be 
written as {Qt = 2/3): 



(10 , as 



.(21) 



TT 



(33) 



The leading order term is 



4'')=4yii7(0,0,xt), 



(34) 

where Xt are given by Eq.([H]) with i = t. At the NLO, assuming an MS top mass we 
have: 



m. 



where: 



Ax 



[C3 - AH{0, -1, 0, x) + H{0, 1, 0, x) - 5/7(1, 0, 0, x)] 

\x(l — x'^'^ 



(35) 



{X - 1)2 

4x [2{1 - x)^ - (2(1 - x^)] 

{X- 1)3(1 + x) 

6x(l + x) 



{x-iy 

6x 



-H{0,0,0,x) - ■ 
-H{0,x) + 



H{0,x) + 



^x(l + X 



{X - 1)3(1 + X) 



H{0,0,x) 



X-1 3 1+x 



Qx 



H{0,0,x). 



(a;-l)(l + x) ' ' ' {x-iy 
The corresponding expression for an OS top mass is given by: 



s^'^'' = si'''^\xt) + ^sr'ixt). 

4.2 Production Cross Section gg ^ A 

The expressions for the relevant quantities in the gg —>■ A production cross section can be 
easily obtained from those in Section [31 with the substitutions: T{Ri) — > 1/2, C{Ri) 
Cp, ^ S, Q ^ K.. In particular, the Born-level partonic cross section (Eq. (l23l) ) is: 
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.(2/,6) , 



(36) 
(37) 

(38) 



(0) 



128 V2 7r 



(39) 



with /C|^'^ = St^''\ The NLO virtual contribution to the gluon fusion subprocess (Eq. fl26|) ) 



is: 



4 ) + ^^''^ 



(40) 



where = (11 Ca — 2nj)/6, with nj the number of active flavor, and 

j^(2i) ^ l^j^my' l^c^/cf + CA/Cf '""^^(xt)^ + h.c. (41) 

with /Cf ''^^^ = and 

^* = (^3Ij^t^3 - Hil, 0, 0, x) - 2n,ix)] + -^—^HiO, 0, 0, x) . (42) 

Eqs.( l38|l4T!) are in agreement with the corresponding expressions presented in Ref. [5]. 

5 Conclusions 

In this paper, we considered the virtual NLO QCD corrections to the processes if — 77 
and gg H. We assumed the couphng of the Higgs boson to the photons and gluons 
to be mediated by fermionic and scalar loops. We provided analytic formulas for these 
corrections that are valid for arbitrary mass of the fermion or scalar particle running in the 
loops and of the Higgs boson. They are given in a very compact form as a combination of 
HPLs. 

The calculation here presented was done using the Laporta algorithm for the reduction 
of the scalar integrals to the Mis and the differential equations method for the evaluation 
of the latters. A part of the Mis needed for the calculation was already known in the 
literature. We explicitly give the analytic results for the Mis that were not known. 

We checked our results for the decay width of the Higgs boson in two photons and for the 
partonic cross section of the gluon fusion by performing an independent calculation in the 
region of small Higgs mass via an asymptotic expansion in the variable r = m\/m? ^ 1, 
with m the mass of the fermion or scalar particle, up to the first 4-5 orders. 

We considered also the NLO virtual QCD corrections to A 77 and gg —>■ A assuming 
the coupling of the pseudoscalar boson to the external particles mediated by a fermion. 

We find complete agreement with the results previously known in the literature con- 
cerning the production and decay of a (pseudo) scalar Higgs boson mediated by fermionic 
loops. This provides an independent check of the formulas given in Refs. [H [5] and extends 
them to the case of a scalar particle running in the loops. 
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Note added 

After our work was completed a paper on a similar subject has appeared on the Web [27] . 
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We did not check yet our formulas against theirs. 



Appendix A: Results in the Low and High Higgs Mass 
Regimes 

In this Appendix we present approximate results that are valid in regions in which the mass 
of the Higgs boson is either much smaller or much larger than that of the particle running 
in the loop. 



A.l i7 ^ 77 

The expansions of Eqs.([5]lZl) "with respect to = mjj/mj, when -C 1 read: 



,(10, 



11 
7 + — r 



19 



30 420 



29 



4200 



41 



34650 



3 90^^ 126^ 



4 7 
3 ~ 90 

12 1 2 
3 ~ 45^^ ~ 140^ 



13 



while in the opposite case, ^ 1 we find: 



12600 
2 



-r 



51975 



1575 



1 



4158 



J-f')(r) = 2+(l2-61n2(-r))^+(241n(-r) + 121n2(-r))-3^ + o(^-3^^ 



(8-21n^(-r)) - - (S In(-r) + 8 In^(-r)) ^ + 



Ti^'\r) = - + 41n2(-r)l + 



(Al) 
(A2) 
(A3) 

(A4) 
(A5) 
(A6) 



The expansions of Eqs. (ll5ll20p when <C 1 are 



'1/2 



■J I 



{2l,b) 



r) 



19 

1 r — 

270 



104 2 6313 



14175 



'60^^ 42^ 



29 



13 

~ 2800' 
4973 



15876000 
o 16 



r" + 



5083 



74844000 



17325 
3137 



-r 



4 216 226800 
1 1 



-r 



4 15 



9 o 2 o 

ny ny ly'--^ 



882000 
5 



560 



525 



5544 



1180367 
2095632000' 



(A7) 
(A8) 



+ 0(r5), (A9) 

(AlO) 
(All) 



while in the opposite regime we have: 



1 



(2l,a). 
/2 I 



36 

36 + :^C'(2) - 4C(3) - 4(1 + C(2) + 4C(3)) In(-r) - 4(1 - C(2)) In^(-r) 
5 



+ ln=^(-r) + — In^(-r) 



144 

48 + 8C(2)-— C'(2)-16C(3) 
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T 



{2l,b) 



+16(1 + 4C(3)) In(-r) + 2(11 - 8C(2)) \Yi^{-r) - 12 In^(-r) 



3 ^ , 



+ 



(A12) 



[12 + 6 In(-r) - 3 hi^^-r)] ^ - [l2 - 24 ln2(-r)] ^ + C (^^^ , (A13) 



1 r 72 

[14 - 3 In(-r)] - + 6 - -('(2) - 24C(3) - 8(1 - C(2) - 4C(3)) In(-r) 
r L 5 



+ (17 - 8C(2)) \n\-r) - | ln3(-r) - I \n\-r) 

6 O 

[6 In(-r) - 3 In^(-r)] ^ + (^) • 



(A14) 
(A15) 



The imaginary part of the expressions in Eqs. (IA4I - IA6P and (1A12I4XT5I) can be easily recov- 
ered via the substitution 

In (— r — ie) ^ In r — ivr . (A16) 



A.2 gg^H 

The expansions of Eqs. fl28ll29p when <^ 1 are 



p{21,Ca) I 
^1/2 



r) 



29 



-r 



3 1080 7560 
1 1 

r + 



1 2 29 

+ 



168000 



6 135 15120 
while in the opposite regime we have: 



1 2 29 

+ 



189000 



Q 3329 A 

74844000 
3 , 1433 4 
29937600 



+ O (r^) , (A17) 
+ O (r^) , (A18) 



12 - ^C'(2) - 12C(3) + 12C(3) In(-r) - (1 + 2C(2)) ln2(-r) 



-— In^(-r) 
12 ^ ' 



1 

- + 

r 



28 + 8C(2) - ^C'(2) - 64C(3) + 8(1 + 6C(3)) In(-r) 



-2(1 + 4C(2)) ln='(-r) - - ln=^(-r) - - In^(-r) 

3 3 



4 + 0^^ 



g(21.C.), 



r) 



6 



+ 



64 



e{2) + 32C(3) - 24C(3) In(-r) + 2(1 + 2C(2)) ln2(-r) 



+ ^ln3(-r) + ^ln^(-r) 
3 6 



(A19) 



(A20) 



The imaginary part of the expressions in Eqs. (lA19IIA20p can be easily recovered via the 
substitution of Eq. (lA16p . 



A. 3 ^4 ^ 77 and gg ^ A 

The expansion of the one- loop result (Eq. fIMP ) when r = m\/mf <C 1 is 



6 45 



1 2 1 

-r — 



280 



1575 



+ O (r^) 



(A21) 
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while for r ^ 1: we have: 



Si^^\r) = 21n2(-r)--81n(-r)4 + C 



The expansions of Eqs. fl36ll37p when r -C 1 read 

,{2l,a),^^ _ 1^ 19 ^2 3533 



151 



-r 



1 

s'^^SeO^ 453600' 141120 

2 



r 



1 1 2 9 o 

ry ry /y'-^ 



4 15 

while in the opposite regime we have: 



560 



525 



-r 



+ O (r^) 



(A22) 

(A23) 
(A24) 



^^C2 - 4C3 + 4(2 - C2 - 4C3) In(-r) - 4(1 - C2) In'(-r) + ln3(-r) 
5 



12 ^ ' 
3 ^ ' 



24 + 8C2 + 32C3 + 8(3 - 2C2) In(-r) - 22 ln2(-r) 



6 ln(— r) — 3 ln^(— r) 



12 - 24 In(-r) 



For Eg. (1421) in the limit r <^ 1, finally, we have: 



t ' 



J- _|_ 



29 



24 60480 



53 



378000 



(A25) 
(A26) 

(A27) 



and for r ^ 1: 



{21,0 a), 



r) 



32 1 

-C2' + I2C3 - I2C3 In(-r) + (1 + 2C2) In'(-r) + - In^(-r) 



1 1 

r 



+ 



28 + 8C2 + 8 In(-r) + 2 ln2(-r) 



J 



(A28) 



The imaginary part of the expressions in Eqs. (1A25IIA26IIA28P can be easily recovered via 
the substitution of Eq. (lA16p . 



Appendix B: Master Integrals for the Two-Loop QCD 
Corrections 

This Appendix is devoted to the analytic expressions of the Mis involved in the calculation. 
They are 11, as it is shown in Fig. [31 

The diagrams (e) and (g)-(k) can be found in j.26j. Diagrams (a)-(d) and (f) are listed 
below. The integrals are performed in Euclidean D-dimensional space and expanded in 
Laurent series of e = (2 — D/2). The mass of the particles running in the loops is denoted 
by m, while /i is the unit of mass of the Dimensional Regularization. The variable x is 
defined in Eq. ([8]). 
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(a) (6) (c) id) (e) (/) 






NWJWwA 



Figure 3: The Master Integrals necessary for the computation of the two-loop QCD correc- 
tions to gg ^ H and H 77. 



A common prefactor of 



2e 



TT 



2(2-e) r2(l + e) 



(Bl) 



is understood and not explicitly shown in the formulas below. 

All the results of this section can be obtained in an electronic form by downloading the 



source files of this manuscript from http://www.arxiv.org 
B.l Topology t = 4 




i=-2 



(B2) 



where: 




-(1) _ 


1 


-2 — 


2' 


= 


1 




2' 


.(1) 


5 


~ 


~2 



X 



{x - ly 



:[4C(3) - 2H{0, 0, 0, x) - 4/7(1, 0, 0, x)] 



2(1 +x) , , 
+ 4 7^H{0,x), 



x-l] 



(B3) 
(B4) 

(B5) 



Fi'^ = -^-2i7(0,l,0,x)+^-^{^C'(2)x+(3a;2-10x + 3)C(3)-2(a:2-l)C(2) 

+[12(x2 - l) + (x - 1)^C(2) -6xC(3)]if(0, x) - [7 - 2x - 9,T2 + 2xC(2)]if(0, 0, x) 
+ [4(a;2 - 1) - 4a;C(2)]/7(l, 0, x) - 12{x^ - 1)H{-1, 0, x) - [3x^ - 8a; + 3] x 



13 



xH{0,0,0,x) + 6{x-IYH{0,-1,0,x) + 2{1 + x'^)H{1,0,0,x) - X 12C(3)/7(1, x) 
+12H{0, 0, -1, 0, x) - 6/7(0, 0, 0, 0, x) - 4/7(0, 0, 1, 0, x) + AH{0, 1, 0, 0, x) 
+24i7(l, 0, -1, 0, x) - 12H{1, 0, 0, 0, x)-8H{l, 0, 1, 0, x) + 8H{l, 1, 0, 0, x)] ] . (B6) 



B.2 Topology t = 5 



where: 




4'^ + O (6) 



(B7) 



^|-^C'(2) - 4C(3)[if(0, x) + 2H{1, x)] + 2if (0, 0, 0, 0, x) 



{x-l) 



+4/7(0, 1, 0, 0, x) + 4i7(l, 0, 0, 0, x) + 8i7(l, 1, 0, 0, a;)| 



(B8) 




(B9) 



where: 



2 17(3) 







X 



C\2) + C(2)i7(0, 0, x) + 2C(2)/7(1, 0, x) + iJ(0, 0, 1, 0, x) 



(x-l)2 

+2i/(0, 1, 0, 0, x) + 3//(l, 0, 0, 0, x)+2if (1, 0, 1, 0, a;)+4i7(l, 1, 0, 0, x) . (BIO) 



For the following MI we have chosen the scalar integral with a scalar product ki ■ k2 
on the numerator, ki and being the two momenta of integration. While the results of 
the other Mis here given do not depend on the routing, F- do. The five denominators 
of the MI under consideration are: Di = [kf + m^], D2 = k\, D3 = [{pi — kiY + m?], 
Di = [{pi - fci + k2f + m% D5 = [ip2 + ki- k2f + m% 



(fci -fca) 



(Bll) 



where: 

F^l = ^, (B12) 
f[^I = 2, (B13) 
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F^'^ = 8 + ^^-l^{3(x2-l)C(2)-[9(a;2-l) + C(2)(l + x2)]i/(0,a;) 



+ U- X + x^)H(0,0,x) - 3x^H(0,0,0,x)} + — ^[4if(-l, 0, x) 

J 2(x — 1) 

+H{1, 0, a;)] - ^^'^^^"'^'^ [2^(1, 0, 0, a;) + H{0, 1, 0, x)] , (B14) 

F^^ = 29 + ^^^1^1^(8x2 -l)C^(2) + 14(x^-l)C(2)- (11 -x + 2x2)C(3) 

-[35(x^ - 1) + (8 - 3a; + a;2)C(2) + (1 - 10x2)C(3)]if(0, x) + [{x^ - 1)C(2) 
+10(1 + x^)CiS)]H{l,x) + [19 + 7x- I8x^ + (4^^ - 3)C(2)]if(0, 0, x) 

-[3{x^-l) - 4(l+x')C(2)]if(l, 0,x) +(l+x2) [c(2)(2iJ(0, -1, x)-H{0, l,x)) 

+H{0, -1, 0, 0, x) + 2/f(0, -1, 1, 0, x) + 2if(0, 1, -1, 0, x) - 2H{0, 1, 1, 0, x) 

+12H{1, 0, -1, 0, x) - 5H{1, 0, 0, 0, x) - 2H{1, 0, 1, 0, a;) + 4iy(l, 1, 0, 0, x) 

+(2a;2-x+8)i/(0,0,0,x) -2(10-5x + 7x2)iy (0,-1,0, x) + 2(l-x + 3x2) ^ 
xH{0, 1, 0, x) - (13 - lOx + 9x^)H{l, 0, 0, x) + 2(9x2 - 2)i7(0, 0, -1, 0, x) 

-9x2iJ(0, 0, 0, 0, x) - (1 + 5x^)H{0, 0, 1, 0, x)| - ^ + [6C(2)g(-l, x) 

J 2[x — 1) 

-34i7(-l, 0, x) + 8H{-1, -1, 0, x) - 3H{-1, 0, 0, x) + 2H{-1, 1, 0, x) 
+2if(l, -1, 0, x) - 2/7(1, 1, 0, x) - 3H{~1, 0, 0, 0, x) - 7H{0, 1, 0, 0, x)] . (B15) 

B.3 Topology t = 6 



J2 + O (e2) 



(B16) 



where: 



.(5) 



4x2 



m4(x - 1)4 
4x2 



-3C(3) - C(2)if(0, x) - 2if(0, -1, 0, x) + HiO, 0, 0, x)] , 
16 

-(2(2) - llC(3)i/(0,x) - 12C(3)if(l,x) - 5C(2)i7(0,0,x) 



(BIT) 



° ~ m4(x-l)H 5 

+4C(2)i7(0, -1, x) - 2C(2)if(0, 1, x) - 4C(2)/7(1, 0, x) + 3i7(0, 0, 0, 0, x) 
+2H{0, 0, 1, 0, x) + 2H{0, 1, 0, 0, x) + 4i7(l, 0, 0, 0, x) + 16/7(0, -1, -1, 0, x) 
-10i7(0, -1, 0, 0, x) - 4if (0, -1, 1, 0, x) - 4:H{0, 1, -1, 0, x) - 8/7(1, 0, -1, 0, x) 



-14i7(0,0,-l,0,x) 



(B18) 
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